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Abstract 

The remarkable strength of glasses is examined using the random first order transi- 
tion theory of the glass transition. The theory predicts that strength depends on elastic 
modulus but also on the configurational energy frozen in when the glass is prepared. 
The stress catalysis of cooperative rearrangements of the type responsible for the su- 
percooled liquid's high viscosity account quantitatively for the measured strength of a 
range of metallic glasses, silica and a polymer glass. 

A fundamental question about solid matter is what ultimately determines its mechanical 
strength. Glasses, in the popular mind, are easy to break but in fact, if surface cracks are 
carefully avoided, glasses turn out to be intrinsically quite strong. Nearly a century ago, 
Frenkel provided an elegant argument for the maximum stress that a solid could withstand 
[]]. Crystalline metals were found to be hundreds to thousands of times weaker than the 
Frenkel estimate [2!]. This observation inspired the extremely fruitful ideas of dislocations 
and grain boundaries that provide easy ways for polycrystalline metals to rearrange and 
plastically deform [3j HI El E] . Glasses come much closer to the Frenkel limit but still fall 
short in strength [7j. In this paper we explore quantitatively the notion that the mechanical 
failure of glassy materials ultimately arises from strain catalyzed rearrangements of the same 
kind as those responsible for the high supercooled liquid viscosity. The idea that there is a 
relation of yield strength to the glass transition itself is not new and has been examined in 
various ways [SI El UHl El [III ETJ]- Here we go further by exploiting the current quantitative 
understanding of cooperatively rearranging regions that has emerged from the random first 
order transition (RFOT) theory of glasses [131 US ESI HSl [T71 HU [19] in order to make 
some specific predictions. RFOT theory describes the microscopic origin of cooperatively 
rearranging regions and predicts they are compact containing a few hundred molecular units 
near the laboratory glass transition temperature T g . These regions become more fractal, 
resembling strings or percolation clusters [20] at higher temperatures where flow is no longer 
thermally activated [21] but rather dominantly collisional. The quantitative predictions 
of RFOT theory concerning the well-established thermodynamic/kinetic correlations in the 
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viscous liquid state, dynamical heterogeneity in supercooled liquids [18J and the aging [22] 
and rejuvenating [T5] properties of the glassy state proper agree quite well with observations 
[25] . It is thus natural to enquire as to what the theory predicts for the material strength of 
glasses. 
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Figure 1: We show in the upper part of the figure schematic snapshots of local rearrangement 
starting from an initial frozen configuration in an imposed stress field a. Following Lubchenko 
and Wolynes [22] the lower left panel shows the spectrum of possible free energy minima for 
a large sample of glass. Levels are listed in order of the internal free energy 0, comprising the 
potential energy along with a vibrational entropic contribution. When the glass is trapped 
in a single such state, local regions of size N can rearrange to new minima while only weakly 
disturbing their environment elastically. Connected energy levels are shown in the next 
two panels. When an imposed stress a is imposed the energy levels are shifted and the 
energy cost of moving N particles is reduced by an amount (na 2 /2G)NVb ea A where G is the 
elastic modulus and k is a factor that includes the elastic response of the environment which 
does not shift to a new minimum. H>ead is the volume of a molecular unit. Eventually for 
sufficiently large N a distinct structure is formed coincident in free energy with the initial 
state, allowing irreversible motions to occur. 



We begin by reviewing how activated events occur in liquids and glasses in the absence 
of stress. The easiest way to conceptualize activated events in the RFOT theory is through 
what is called the landscape library construction by Lubchenko and Wolynes [22]. This 
construction has also been used to define point-to-set correlation lengths J2U [25] allowing 
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Figure 2: The local level libraries in Fig. [T] when thermally averaged yield a free energy 
profile for rearranging a specific region as a function of the number of displaced particles 
N. The average mismatch energy is balanced against a term containing the configurational 
entropy (from averaging over all the states), any initially excess energy frozen into the glass 
along with a contribution from relaxing strains via reconfiguration in an imposed stress field. 
The activation barrier is lowered by the imposed stress, eventually vanishing when the stress 
is sufficiently large, leading to rapid failure of the glass's structural integrity. 

many key points of RFOT theory to be confirmed via computer simulations [26j [27J [28j [29] . 
This construction is schematically pictured in Fig. [TJ In mean field RFOT theory, below a 
dynamical transition temperature T4, the system becomes trapped in one of an exponentially 
large number of possible metastable states which are minima of a free energy functional [14J. 
For molecular fluids these states can be taken as nearly structurally synonymous with the 
inherent structures which precisely correspond to minima of the potential energy [16J, but the 
individual stability of these states at finite temperature depends not only on their energy but 
also on their vibrational entropy. Irreversible reconfiguration events eventually take place by 
rearranging molecules in ever larger regions of size N until a critical size is reached. Above 
the Kauzmann temperature, Tg, the configurational entropy is extensive and so as the size 
of a reconfiguring region increases the number of possible local rearrangements grows as 
well. Generally moving to any one of these rearranged structures costs free energy because 
the environment of the rearranging region does not fit the new locally accessible alternative 
structures as well as it fits the original free energy minimum from which rearrangement 
starts. The typical mismatch energy AE(N) near the Kauzmann transition scales as ^N x . 
The power law in mean field theory represents a surface energy [T3] so the exponent x = 2/3 
but scaling arguments [13] suggest there should be a somewhat weaker scaling with x — 1/2 
near an ideal glass transition at TV due to wetting from the numerous alternative states that 
can be interpolated between the fixed environment and the core of the rearranging region. 

Xia and Wolynes showed the coefficient in the mismatch energy can be computed near 
TV by assuming a locally sharp interface and by making a microscopic estimate using density 
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functional theory of the localization free energy which is entropic: 

Sv 7 ^, „ . / a L a 2 \ 

Here oil determines the size of the vibrational fluctuations in a minimum and is roughly 100, 
reflecting displacements following Lindemann's stability criterion allowing localized motion 
of about one tenth of the interparticle spacing, a . 

Above Tk any mismatch energy can however be overcome by the entropic driving force 
favoring reconfiguration to one of the many alternate structures, -Fbuik(^) = —Ts c N where 
s c is the configurational entropy per particle. Balancing i*buik(-W") an d AE(N) gives an ac- 
tivation free energy to be overcome for irreversible rearrangement, AF$, which is a function 
of s c . AF* diverges near Tk as s c vanishes. This prediction then connects the kinetics of 
rearrangements with thermodynamics, a hallmark of the RFOT theory. Using the approxi- 
mate coefficient 7 obtained by Xia and Wolynes the absolute magnitude of barriers is also 
predicted to follow an Adam-Gibbs-like relation AF^ = A/s c but with a specific numerical 
value for A = ^A; B (ln(a L a 2 /7re)) 2 . Because the Lindemann parameter q £ depends only 
weakly on the potential, in RFOT theory then AF'/^T is again dominantly a function of 
the configurational entropy, across a range of substances. 

The landscape library argument can also be used in the so-called "aging" regime to de- 
scribe motion in the glass [22]. In the aging regime, the initial configuration is not one chosen 
from the thermal equilibrium ensemble at the ambient temperature but instead structurally 
resembles a system that was equilibrated at a higher so-called "Active" temperature. In a 
simple quench to low temperature the Active temperature initially is the laboratory glass 
transition temperature T g . For this nonequilibrium situation the initial configuration then 
will not only gain entropy by reconfiguring locally but also will release an additional energy 
per particle A$ which represents the energy frozen in at the glass transition [22]. If we 
assume the configurational heat capacity has the empirical form Ac p (T g ) ■ (T g /T) this excess 
energy is A$ = Ac p {T g )T g In ^. 

Owing to this excess driving force, reconfiguration events occur sooner in a glass than 
they would in a liquid structurally equilibrated to the same ambient temperature. The 
nonequilibrium activation free energy can be written in terms of the activation free energy 
for an equilibrium liquid at a higher specific configurational entropy: 

AFi e =AFi q (s c +^) (2) 

AF£ Js c ) is the function giving the activation barrier in the liquid, written in terms of 
its configurational entropy. In the glass below its T g this formula implies the rates follow 
something close to an Arrhenius law but with an activation enthalpy diminished from what 
it was at the higher temperature at which it fell out of equilibrium. In this way, this RFOT 
argument accurately predicts the so-called nonlinearity parameter x describing the ratio of 
activation enthalpy for motion in an equilibrated liquid to that for glasses that have fallen 
out of equilibrium [22J. 
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The Volger-Fulcher law, while describing the deep glassy behavior, breaks down at higher 
temperatures in supercooled liquids. In mean field theory this breakdown occurs because the 
mismatch coefficient 7 itself vanishes at the mean field spinodal temperature Ta [EH EHl ED] • 
Schmalian, Stevenson and Wolynes have argued that the Volger-Fulcher relation will actually 
break down at a somewhat lower temperature T c because the shape of correlated activated 
regions changes at higher temperatures in such a way that the mismatch energy now scales 
linearly in N [20] • The high temperature, entropically favored shapes are lattice animals 
whose exposed surface scales directly with their number of constituents as does their shape 
entropy. In this regime the scaling of mismatch energy can be written as Vi nt b where b is 
the number of equivalent broken bonds at the surface of the rearranging region. Schmalian, 
Stevenson and Wolynes showed that near Tk, b is approximately 3.2N and the coefficient of 
Vi n t can (like the surface mismatch energy 7) be obtained from density functional reasoning 
v int = (l/z)(3/2)k B T\n(a L a 2 /ne). 

The free energy profile for such fractal rearranging regions either monotonically increases 
with N or decreases monotonically with N. This means that there will be a change in the 
rearrangement mechanism from activated dynamics to one dominated by collisions at high 
temperature. The crossover to barrierless reconfiguration is thus determined by the condition 
s c (T c ) = sj? erc — A$/T c = Ac p (Tg)Tg/T K (l — T K /T C ) . For equilibrated liquids this relation 
predicts crossover temperatures agreeing with those found using Stickel plots [31]. The 
RFOT argument then also predicts barrierless reconfiguration will occur for nonequilibruim 
glasses if heated when 

k B T c k B T K ( s^ c T K \- x 

Ac p (Tg)T g - A$ Ac p (T g ) T g \ Ac p Tg ) 

This specific prediction for crossover to collisional dynamics in superheated glasses has not 
yet been tested in the laboratory although it would be interesting to check it experimentally 
by using lasers to rapidly heat glasses. 

What does RFOT theory then suggest about how reconfiguration events occur under 
stress? If a sample of glass is put under a uniform shear stress, a, the energy per unit volume 
is immediately raised by an amount a 2 /2G where G is the elastic shear modulus [32]. It has 
long been known this energy can be explosively released by cracking the glass. This effect is 
demonstrated by the famous Prince Rupert's drops [33]. The stress need not always lead to 
cracking directly. It is reasonable to expect that, of the myriad of possible states envisioned 
in RFOT theory, a significant fraction will also allow this stress energy to be released without 
cracking or forming voids. Indeed a vanishing stress energy of the rearranged state is expected 
since the most stable mean field free energy minimum corresponding to delocalized molecules 
can be thought of as being a typical disordered liquid ensemble and is thus completely 
incapable of sustaining static shear. Since the imposed stress energy can be removed by 
appropriately rearranging a region, imposed strains will lower the activation barrier and will 
catalyze the rearrangement. If the stress is sufficiently large the rearrangement may even 
occur without any significant barrier at all, just as takes place at the thermal crossover at 
T c . This crossover to barrierless reconfiguration would thus give the limiting strength of 
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the glass if we assume there are no easier routes for the glass to rearrange (like cracking). 
Strain catalysis means that a stressed glass will always flow at some finite rate even for 
the smallest stresses [HI E2] and thus a glass will deform, given time, at somewhat lower 
stresses than this limit. This gentler flowing situation is probably quite relevant in many 
practical situations. Flow itself can act to further catalyze rearrangements. The resulting 
additional enhancement of reconfiguration speed is a facilitation or mode coupling effect. 
Lubchenko has shown that this effect (that would be contained in a more complete RFOT 
theory including mode coupling effects [31]) does a good job describing the crossover to 
steady state nonlinear rheology [35]. Similar effects have been studied in mode coupling 
treatments of dense colloid rheology [361 EZ]- We will, in this paper, however, concentrate on 
the immediate effect of stress on the activated events that occur before flow starts and leave 
the physics of developed plastic flows for future work. The limiting strengths we calculate 
in this paper then should be upper bounds representing extremely rapid failure of the glass. 





Figure 3: The upper panel shows a uniform strain field acting on the glass sample in its 
original state. A fluidized region allows the surrounding material to elastically deform in a 
nonuniform way in the imposed stress field shown schematically in the lower panel. This al- 
lows additional strain energy to be released without costing any additional mismatch energy. 

Naively speaking, in order to compute the effect of stress catalysis on the activation bar- 
rier one merely needs to account for the strain energy lost in the fluidized region and thus 
must add to the bulk thermodynamic driving force term (— Ts c — A$)iV an additional con- 
tribution a 2 NVbead/2G to compute the lowering of the thermal barrier for compact clusters 
or to find the limiting stress where barrierless rearrangement may occur. Here Vb ea d is the 
volume of a "bead" i.e. movable unit of the glass, which can be inferred from the molar 
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fusion entropy [301 123]. There is a subtlety, however; as pictured in Fig. |3j fluidizing a 
region of the glass also allows Hookean elastic rearrangements of the surrounding matrix to 
occur without it being necessary for the outside region to move to any alternate free energy 
minima. This outer region while elastically responding thus does not elicit any mismatch 
energy. The strain energy relieved by rearrangement of a region of size N (38] nevertheless 
becomes larger than a 2 NVb ea d/2G. 

For an arbitrarily shaped rearranging region the exact calculation of the additional strain 
energy relieved by harmonically distorting the outer region would seem to be a complex 
problem in elastic theory. The result for spherical regions, however, has been known for some 
time where it has been used to develop the theory of the elastic modulus of composite media 
containing holes [39]. The calculation is mathematically quite analogous to the calculation 
of intrinsic viscosity made first by Einstein [10] for the effect on viscosity of suspending solid 
colloidal particles in a fluid and still more closely follows Taylor's analysis of the viscosity 
changes due to suspending liquid drops or bubbles in a fluid |4l|. For spheres the additional 
energy released (analogous to the excess viscous dissipation in the hydrodynamic problems) 
is still proportional to the sphere's size and according to MacKenzie depends on the Poisson's 
ratio characterizing bulk versus pure shear deformations. Taking over MacKenzie's correction 

2 

gives then an energy increment for rearranging a region of size N, Ai? e i astic = ftf^iVVbead 
where k = 3 — 6/(7— 5u) in terms of the Poisson's ratio v. For the typical Poisson's ratio of 
metallic glasses k ~ 1.8. 

Directly calculating the correction for non-spherical shapes is indeed complex. In addi- 
tion, considerable distortions from spherical geometry are energetically still more favorable 
in relieving stress than for the region to remain compact and thus at high stress such distor- 
tion should again lead to barrierless breakup, just as critical flow rates lead to dissolution 
of drops in emulsions [121 H31 HI]. The latter problem has led to an extensive literature 
[151 [38] . We will use, nevertheless, the spherical value of the correction k for all shapes of 
cooperatively rearranging regions. We suspect this simplification is probably not too bad 
for small stresses and not too far from T K . This surmise is buttressed by the experience for 
the corresponding hydrodynamic problem of computing the intrinsic viscosity of complex 
shapes, a problem that has been extensively studied in polymer chemistry [38]. In that 
problem the shape effects are quite modest. By adding the increased relieved strain energy 
to the reconfiguration driving force, in analogy with equation ([2]), the activation barrier for 
flow in a strained glass can again be written in terms of the function giving the barrier 
for equilibrated liquids AF t = AF t (s c + A$/T + Ka 2 V bead /2GT). With this simplification 
then barrierless reconfiguration should finally occur when s c (T c ) = sj? erc — ^ — K^fVbead ■ 
As in the popular J-point scenario [16] an apparent spinodal to reconfiguration apparently 
can be approached by tuning either the temperature T or the stress a. 

The argument relating barriers in the glass under stress to those for thermal motions in 
the equilibrated liquid should hold for temperatures not too far from Tk since the shapes of 
rearranging regions are determined entropically. There are corrections, however, away from 
Tk ■ At very high temperatures near to the mean field spinodal T4 the mismatch energy cost 
goes down, leading to an additional weakening of the glass. Conversely at low temperatures, 
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much below Tr-, we also must account for both the fact that the mismatch energy becomes 
pinned at its Tk value and that at the same time the importance of shape entropy is lessened 
by the diminished temperature. A detailed account of the latter effects is contained in the 
Supplementary Information. When the latter effects are included along with the calculation 
of the excess energy we find an explicit equation for the limiting strength a*: 

1 2Gk B T 77 T~ K 1 Ac„(T ) T Q ' T g \ 

The contribution in this expression involving Ac p represents the weakening caused by the 
excess energy which has been frozen in at the glass transition. If we could be cosmologically 
patient this excess energy would disappear by annealing to T g = Tk giving then the ultimate 
achievable strength of a glass. At very low temperatures the strength of this most stable glass 
will be then a* deal = ^3.QGksTK jVbead ■ We can write the shear modulus in terms of the 
spinodal temperature Ta and the bead size, following Lubchenko by estimating vibrational 
displacements from Debye continuum theory and assuming a limiting Lindemann ratio for the 
maximal thermal excursions (471 Using this relation gives G = 24.9fcsTA/Head- If we now 
take Ta ~ T m the melting point and use the typical ratio of T^/T m of between 0.4 and 0.7 we 
find the ultimate limiting a* is proportional to the elastic modulus. Such a linear correlation 
between strength and modulus, resembles Frenkel's estimate, and indeed has been examined 
experimentally. We find the ideal limit strength from RFOT theory to be uniformly about 
30% higher than Frenkel's. The weakening of the glass due to energy frozen in at the glass 
transition is however substantial. This excess energy lowers the strength quite a bit below 
the RFOT ideal value and below the Frenkel value but still gives strengths greatly exceeding 
the measured strength of polycrystalline metals. We have gathered from the literature data 
for the input thermodynamics. We then compared the RFOT theory predictions to measured 
strengths for some metallic glasses, silica and a polymer glass, PMMA. Details of the input 
data and measured strength data can be found in the Supplementary Information. 

In Fig. Ha] we display results for the strength versus shear modulus. The predicted 
strengths generally exceed but are close to the measured yield strengths. On this plot we 
also show both the Frenkel estimate and our Tg = Tk ideal value. A typical polycrystalline 
material value of one-hundredth of the Frenkel value is also plotted. Clearly the present 
RFOT predictions account very well not only for the trends but even the actual magnitude 
of the strength. In Fig. [4b] we show the comparison of measured strengths against the 
complete predictions. Not only the glassy metals but also silica and the plastic PMMA have 
strengths not terribly far off the RFOT predictions. Although the main dependence on elastic 
modulus is clear, the RFOT theory results also depend on other quantities, such as Ac p and 
the ratio of the ambient temperature to glass transition temperature. However, as we can 
see, both the predicted ratios of strength to modulus for the measured systems and the 
measured ratios show no overall trend with liquid fragility or glass transition temperature, 
see Fig. SI in the Supplementary Information. Of course since modulus and T g are well 
correlated the absolute strengths themselves do correlate with T g . It may be possible to test 
the theory further. Again, rapid heating should lower the yield strength in a predictable 
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(a) Strength as a function of shear moduli (b) A plot of the strength a*, measured from the 

experiments versus theoretical estimation 



Figure 4: The experimental strength (red triangle) and the predicted strength (black circle) 
have nearly the same slope and are quite different from Frenkel strength (blue solid line). 
Typical value of crystal strength (green dashed line) and strength in the limit T — > 0,T g — > 
Tk (black dashed-dotted line) are also shown in comparison. 



way. In addition, superstable glasses can be made via vapor deposition |48j. Their effective 
temperature corresponds to being roughly half way to the Kauzmann temperature. We see 
their strength should thus be proportionately closer to the Frenkel limit. 

RFOT theory accounts well for the measured strength of laboratory glasses of various 
composition. The good agreement between theory and experiment suggests that the corre- 
lated rearranging regions responsible for high temperature viscosity in supercooled liquids 
also limit the strength of nonequilibrium glasses. There seems to be no necessity to invoke 
then any additional defects of a point-like or line-like character to play a prominent role in 
weakening glasses that are prepared in an ordinary fashion by cooling a melt. 
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1 Derivation of equation (4) 

The free energy cost for cooperatively rearranging a region (CRR) of N sites with b boundary 
interaction is 

F(N, b, a) = [f eq (T) - <f) m (T, <t)]N + v mt b - k B T In Q(N, b) (SI) 

where Q(N, b) is the number of lattice animals of given N and b, and s c is the configurational 
entropy per site. The initial nonequilibrated free energy has an additional elastic energy due 
to stress 

2 

</> in (T,a)=<fi n (T)+K^V bead (S2) 

where (f)° in is the sum of the energetic and vibrational entropic contribution to the initial 
nonequilibrated state. 

The equilibrated free energy at temperature T above T K is 

f eq {T) = <Pk- f dT's c {T') (S3) 
Jt k 

= ^-A Cp (r ff )rJ^^-ln^ (S4) 

where we have used Angell's empirical form: Ac p (T) = Ac p (T g )(T g /T), and s c (T) = 
dT , ACp ^P = Ac p (T g )T g (^L-ij. Note that f eq (T g ) = <f> g - T g s c {T g ), the ideal 
glass state energy is equal to 

( f )K = ( f )g -Ac p (T g )T g \n^. (S5) 
14 



Consider the first term on the right-hand side of equation (1ST]) and let $ n (T) be the bulk 
energy at T g . Then use the relation in equation (1ST)]) to obtain 



feq(T) - <f>in(T,cr) 



-AC p (Tg)Tg 



Tk T 



K ~2G ^ bead 



(S6) 



Substitute this result back into equation (JST 



F{N,b,a) 



-AC p (Tg)Tg 



Tk T 



a 



N + v mt b-k B T\nQ(N,b).(S7) 



If T is below Tk, the excess energy is frozen at the state T = Tk and the configurational 
entropy vanishes. 



F(N,b,a) 



-Ac p (T g )T g In 



T 



a 



— K- 



T K 2G 
In percolation clusters p], for large N, 

(a + ir +ix x 



V, 



bead 



N + v int b-k B T\nQ(N, b). 



(S? 



perc 



cxp 



N 2 <t> 



(S9) 



where a = t/N, and t is the number of unoccupied sites bounding the occupied cluster. Take 
the exponent <p a ^ mean field value of 1/2 and a lattice constant B = 1.124. The number of 
bonds b is directly related to t and should linearly depend on coordination number, z, b ~ 



1.68te/zgc[2]. Recall that the surface energy v. 



hit 



-T K I -kn In 



B 



ire 



3.6907- 



B-LK 



The free energy in equation ( 1S8I) is now a function of A" and t. Minimize this function with 
respect to t, the most probable value of t is i = aN, where a = 3.10 at T = T#. At this most 
probable value, f2 perc becomes a simple exponential function, Q peTC ~ A^, where A = 6.82. 
Each term in equation flS8j) is now proportional to A" and we can write 



F(N,b,a) 



T a 2 

-Ac p (Tg)T g ln JT~ K^QV bead 

+v int l.m^aN - k B TN In A 

zsc 



N 



(S10) 



k B TN 



ACp(Tg) Tg 



k 



B 



T 



T 
ln^ 



1 /ta 2 
fc B T 2G 



6ea<i 



1.68^^a-lnA 



k B T ' z sc 

At the thresholding stress a* where F(A', a*) = 0, one finds 



(sn; 



a 



'2Gk B T 

KVh 



bead 



k B T 



zsc 



'- a — In A 



Ac p (Tg) Ti ln Tg_ 
T T K 



k 



B 



(S12) 
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or 




2Gk B T v int 

KV bead k B T 



1.68 



Zf.cc. 



zsc 



a — In A + 



A Cp(Tg) Tg Tg 1 k B T 



k B T T K Vi n t 



(S13) 



Substituting the numbers in the previous paragraph, we finally obtain equation (4) 



2 Experimental data and numerical prediction 

In Table SI we summarize the input thermodynamic data and measured strengths as well 
as their sources in the literature. The bead count is obtained as described by Lubchenko 
and Wolynes [3] and Stevenson and Wolynes [2] using the melting entropy. All strength 
measurements were all made at room temperature 300 K. 




(S14) 
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Glasses 


P 


AH M 


Nbead 


T K 


T g 


Tm 


Ac p 


G eX pt 


°cxpt 


n* 
prcd 


Refs. 


PMMA 


1188 


4.64 


0.84 


337 


372 


397 


6.68 


3.10 


0.07 


0.08 


H El El El El El m\ 


Si0 2 


2648 


9.6 


0.78 


876 


1500 


1995 


1.67 


31.40 


1.1 


0.90 


[m H2i H31 in us] 


M g 6 5 Cu 25 Y 10 


3978 


8.65 


0.66 


325 


424 


771 


2.41 


19.60 


0.8 


1.29 


[ISlH71II51[igil2Dll2IlB21l25] 


Mg 65 Cu 20 Zn 5 Y 10 


3284 


7.77 


0.79 


325 


404 


702 


2.38 


23.60 


0.88 


0.92 


PS 123 


M g 8 Cu 10 Y 10 


- 


7.21 


0.69 


- 


427 


746 


2.73 


- 


0.8 


0.96 


[mess ill] 


La 55 A1 25 Ni 20 


6140 


7.48 


0.75 


337 


491 


712 


2.15 


13.40 


0.6 


0.52 


[2T1 1271 1251 H51I22] 


La 55 A1 25 Cu 5 Ni 15 


6050 


7.51 


0.82 


328 


472 


660 


2.21 


- 


0.6 


0.47 


[2DCEBJEI] 


La 55 Al 25 Cu 10 Ni 10 


5930 


6.84 


0.74 


332 


467 


662 


2.40 


- 


0.6 


0.48 


[2THT81I2T] 


La 55 A1 25 Cu 15 Ni 5 


6370 


7.21 


0.78 


320 


459 


663 


2.14 


- 


0.6 


0.52 


[271IT81I2T] 


L a 5 5 A1 25 C °5 Cu 10 Ni 5 


6000 


6.09 


0.66 


363 


466 


661 


2.84 


15.60 


0.8 


0.67 


[271 US EU E9] 


Pd 40 Ni 40 P 20 


8951 


7.39 


0.60 


487 


570 


884 


3.28 


39 


2.19 


2.33 


13111 [T8l I3T1 [2T1 122] 


P 40 Cu 30 Ni 10 P 20 


9300 


6.84 


0.61 


497 


593 


798 


3.39 


33.00 


1.8 


2.30 


PHI E2J CLEl EH ESI E2] 


Pd 77 Cu 6 Si 17 


10400 


8.55 


0.58 


553 


637 


1058 


2.74 


31.80 


1.5 


3.11 


[331 EH da [35] 


Zr n Cu 27 Ni 8 Ti 34 


6850 


11.3 


0.72 


537 


671 


1128 


2.25 




2.2 


3.20 


P31 EH E3 EH ES] 


Zr 41.2 Ti 13.8 Ni 10 Cu 12.5 Be 22.5 


6125 


8.2 


0.63 


553 


620 


937 


4.03 


34.10 


2.12 


2.26 


[ISl EH US EH ED SB ES] 


Z r 46.25Ti8.25^ioCu 7 5 Be 27 5 


6014 


9.4 


0.57 


550 


622 


1185 


3.60 


34.70 


1.83 


2.68 


dsi nzi imi esi na sd Ha esi 


Zr 52.5 Cu 17.9 Ni 14.6 A1 10 Ti 5 


6730 


8.2 


0.55 


633 


675 


1072 


4.35 


32.12 


1.88 


2.48 


[2H EH UM EH HOI S3] 


Zr 57 Cu 15.4 Ni 12.6 A1 10 Nb 5 


6690 


9.4 


0.62 


656 


682 


1091 


3.26 


30.80 


1.8 


2.69 


[2D E3 US m H3] 


Zr 65 A l7.5Cu 275 


6744 


12.8 


0.80 


517 


666 


1150 


2.24 


33.00 


1.7 


2.16 


[2DE3EH] 


Table SI: Experimental dat 


a and 


theoretical results of 19 g\ 


lasses: 


density p (kg/ 


m 3 ), 


latent 


heat of fusion AHm 



(kJ/mol K), number of bead Nb ea d, Kauzmann temperature T% (K), glass transition temperature T g (K), melting tem- 
perature Tm (K), heat capacity change at glass transition temperature Ac p (J/mol K), experimental shear modulus G expt 
(GPa), measured strength a* t (GPa), and theoretical estimated strength a* d (GPa). 
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Figure SI: The ratio between strength and shear moduli versus heat capacity change at T g 
and the glass transition temperature T g . The black circles are the RFOT theory predictions 
and the red circles are the measured values. 



We also show here in Fig. [S2] the comparison of measured elastic moduli with those 
predicted via the relation G ca i = 2A.Qk B T m /V^ e& ^ from thermodynamic data along with the 
Lindemann relation and the assignment Ta ~ T m . 
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Figure S2: Comparison between the measured elastic moduli and those calculated using the 
Lindemann criterion. 
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